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Reaction theory is an essential ingredient when performing studies of nuclei far from stability. 
One approach for the calculation of breakup reactions of exotic nuclei into two fragments is to 
consider inelastic excitations into the single particle continuum of the projectile. Alternatively one 
can also consider the transfer to the continuum of a system composed of the light fragment and the 
target. In this work we make a comparative study of the two approaches, underline the different 
inputs, and identify the advantages and disadvantages of each approach. Our test cases consist of 
the breakup of ^’^Be on a proton target at intermediate energies, and the breakup of ®B on ®®Ni at 
energies around the Coulomb barrier. We find that, in practice the results obtained in both schemes 
are in semiquantitative agreement. We suggest a simple condition that can select between the two 
approaches. 

PACS numbers: 24.10.Ht, 24.10.Eq, 25.55.Hp 

I. INTRODUCTION 

A large fraction of present nuclear physics encompasses 
the study of nuclear structure far from stability using re¬ 
action measurements. Extraction of fundamental struc¬ 
ture therefore requires adequate reaction theories. As 
one of the most important tools, breakup offers a unique 
opportunity to benchmark reaction theories. From the 
early days it became clear that the standard models to 
breakup of stable nuclei needed revision Q. Many of 
the lessons learnt from the deuteron breakup have since 
become a source of inspiration for the rare isotope re¬ 
action community 0. Different groups designed various 
reaction models, tailored to specific systems and using 
particular approximations. Albeit the variety, the state 
of the art of the existing reaction-model panorama has 
become increasingly unsatisfying: at present we already 
have a handful of models that produce results for a spe¬ 
cific case but we are missing a general effort of a consis¬ 
tent comparison between the various approaches. In the 
few cases where two different models are applied to the 
same problenn there is often a disparity in the predictions 
1^ ■ It is timely to make the necessary links 
between the available models. Some studies, comparing 
approximations such as eikonal, adiabatic, local momen¬ 
tum, have been recently performed Here, we 

work within a framework where no such approximations 
are present. 

Roughly speaking, current breakup reaction theories 
can be divided into two main categories. On one side, 
some methods model the breakup process as an excita- 
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tion of the projectile to the continuum spectrum of the 
projectile. This is the case of the CDCC (Bin ap¬ 
proach. On the other side, some other methods, such 
as the semiclassical transfer to the continuum developed 
by Brink and Bonaccorso or the 

post-form DWBA approach [l^ 1^ 1^ 1^ . treat the 
breakup process as the transfer of one of the fragments 
to the unbound states of the target. Intuitively, it is 
obvious that both continua do correspond to the same 
three-body continuum, expressed in different coordinates 
systems. However, it is not clear to what extent this 
equivalence is fulfilled in a practical calculation. In this 
work we try to shed some light on this problem by ap¬ 
plying both approaches to the same reaction and using, 
whenever possible, the same physical ingredients. 

Although many kinds of observables could be calcu¬ 
lated in both methods, we make special emphasis on core 
energy and angular distributions, since these observables 
are particularly important in currently measured reac¬ 
tions with radioactive beams. 

The paper is organized as follows. In section H we 
briefly review the three-body breakup and transfer to the 
continuum approaches, in the form used in this work. 
In Sec. HI, we apply these formalisms to the reactions 
p+^^Be and ®B+®®Ni. A discussion of these results is 
presented in Sec. IV and conclusions are drawn in Sec. 
V. 



(1) (2) (3) 

Figure 1: (Color online) The three Faddeev components for 
the problem of a two-body projectile (c -|- x) impinging on a 
target T. 
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II. THREE-BODY REACTION MODELS 

When considering reactions with light radioactive 
beams, it is customary to model the incoming projectile 
as a two-body system (in fact sometimes the projectile 
has a clear three-body structure and models that handle 
three-body projectiles are underway). In principle, the 
solution to the reaction problem can be obtained exactly 
from solving the Faddeev equations with the appropriate 
boundary conditions. In the Faddeev formalism |2.*ll^ . 
the three-body wavefunction is written as a sum of three 
Jacobi components represented in Fig. ^ Each compo¬ 
nent is defined by the intercluster coordinate between two 
of the subsystems ({r^; z = 1,2,3}) and the relative coor¬ 
dinate of this pair to the third cluster ({Rz; z = 1, 2, 3}). 
Asymptotically, the Faddeev component i contains con¬ 
tribution from the bound states associated to the pair 
with relative coordinate r^, plus a contribution coming 
from three-body breakup. Therefore, while rearrange¬ 
ment channels are confined to specific Faddeev compo¬ 
nents, breakup is distributed among the three compo¬ 
nents. Consequently, extracting the breakup observables 
requires complicated transformations among Jacobi coor¬ 
dinates. Besides, solving the Faddeev coupled equations 
is a very difficult task, specially for three charged parti¬ 
cles where Coulomb plays an important role. Inclusion 
of absorption in the intercluster potentials, which is re¬ 
quired when the subsystems have internal degrees of free¬ 
dom is also an open problem, although some promising 
work is already in progress 

It is well known that basis states belonging to different 
Jacobi sets are not mutually orthogonal. Furthermore, 
for each Jacobi set, a complete basis of the three-body 
bound and unbound spectrum can be constructed. Then, 
it could be possible, in principle, to describe the reac¬ 
tion observables of a thee-body scattering problem us¬ 
ing uniquely states from one of the Jacobi components. 
This is in fact the procedure followed by the Continuum 
Discretized Coupled Channels (CDCC) formalism. This 
method has been applied for more than two decades to 
the scattering of weakly bound (two-body) projectiles by 
light and heavy targets. For the scattering of a compos¬ 
ite projectile A{=c+x) by a target T, the CDCC method 
defines the model three-body Hamiltonian: 

H = Krel + Hint + Ux + Uc, 

Hint Hint + Vxc ( 1 ) 

where Krei is the kinetic energy for the projectile-target 
relative motion, Kint is the internal kinetic energy of the 
projectile, Ux and Uc are the x — T and c — T interac¬ 
tions and Vxc is the x — c binding potential. As to the 
internal degrees of freedom of the target, in the standard 
CDCC method, only the target ground state is considered 
explicitly. Therefore, the fragments are not allowed to 
engage in arbitrary processes with the target. For exam¬ 
ple, processes in which one of the dissociated fragments 
is absorbed by the target, or in which the target inter¬ 
nal degrees of freedom are excited, are excluded from the 


model space. Also, rearrangement channels correspond¬ 
ing to cluster-target bound states are by construction 
excluded from the CDCC model space and hence, those 
observables associated with these two-body channels can 
not be obtained from the asymptotics of the CDCC three- 
body wavefunction. The model space spanned by CDCC 
allows only the calculation of elastic breakup and leaves 
out those processes related to inelastic breakup. 

In order to take into account the effect of the excluded 
channels, the interactions Ux and Uc are usually taken as 
phenomenological optical potentials obtained, for exam¬ 
ple, from the fit of the elastic data at the same energy 
per nucleon. By contrast, the interaction Vxc is taken to 
be real, and chosen to reproduce known bound and/or 
excited states separation energies, or resonance energies. 

The full three-body space is truncated by setting a 
maximum excitation energy for the projectile. Moreover, 
the c — X relative angular momentum is also restricted 
by considering only a limited number of partial waves. 
In order to deal with a finite set of coupled equations, a 
discretization of the continuum states into energy inter¬ 
vals (bins) is also performed. This procedure should be 
regarded as a practical method of making the problem 
numerically solvable, rather than an additional approx¬ 
imation. In fact, it has been shown that the calculated 
observables are essentially independent of the method of 
discretization |2^ . 

Within this restricted model space, the three-body 
wavefunction is expanded in eigenstates of the internal 
Hamiltonian Hint as 

Ni 

^'W(ri,Ri) = ^(/)«(ri)xa(Ri), (2) 

a—0 

where Ni is the number of states considered, a represents 
all angular momentum quantum numbers as well as exci¬ 
tation energies of the projectile, </a(i") are the eigenstates 
of the two-body Hamiltonian Hint and Xq(R) describes 
the relative motion between the projectile A = c + x and 
the target T. This expansion of the three-body wavefunc¬ 
tion is inserted into the Shrodinger equation that, when 
projected into the considered internal states, provides a 
set of coupled equations. 

Within the CDCC scheme, the breakup process is 
treated as inelastic excitations of the projectile A into 
the continuum c + x due to the interactions with the tar¬ 
get T [11 0 m. A pictorial representation of these 
couplings is given in Fig.[2Da). The couplings responsible 
for this excitation, as well as the diagonal potentials, are 
obtained by folding the phenomenological interactions Ux 
and Uc with the internal wavefunctions, i.e. 

?7q,;q,/(Ri) = {4>a\Ux + Uclff’a') ■ (3) 

Applications to the breakup of ®B at low and inter¬ 
mediate energy regimes have been very successful in de¬ 
scribing the data P.13.12^. 

Unlike the Faddeev method, the CDCC approach uses 
only one of the three possible sets of Jacobi coordinates. 
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As noted above, rearrangement channels corresponding 
to cluster-target bound states are not part of the CDCC 
model space and, therefore, the CDCC three-body wave- 
function is not adequate to predict observables associ¬ 
ated with these two-body channels. On the contrary, it 
has been argued that, provided that the model space is 
sufficiently large , the total three-body breakup is 

contained in the CDCC wavefunction and, hence, can be 
extracted from its asymptotics. 

Although considerably simpler than its Faddeev coun¬ 
terpart, solving the CDCC problem is also a complicated 
task. In some cases, particularly with heavy targets, long 
range interactions usually lead to convergence problems. 
An additional difficulty is that, in many breakup exper¬ 
iments, scattering observables (differential energy cross 
sections, angular distributions, etc) are obtained with re¬ 
spect to one of the projectile fragments (this is indeed 
always the case of inclusive reactions). Given the choice 
of coordinates, the CDCC observables are more naturally 
expressed in terms of the projectile center of mass, and 
its internal excitation energy. Converting to one of the 
fragment’s coordinates requires a complicated kinematic 
transformation j^. Furthermore, due to the restricted 
model space, the CDCC description is not expected to 
be good in the region where channels outside the model 
space play an important role. Discrepancies in large an¬ 
gle scattering data, observed in early applications of the 
CDCC method to deuteron and ^He breakup [^, have 
been attributed to this fact. 

A way to circumvent the two latter criticisms is to use 
the T-matrix formalism in post-form and approximate 
the incoming exact three-body wavefunction appearing 
in the exact scattering amplitude: 

Tpost = (xisVirV.c + U,- (4) 

by the CDCC wavefunction, i.e., 4'-''’^ w -q/CDCC_ 

equation, B = x+T, xi's is the distorted wave generated 
by the (arbitrary) distorting potential C//(R') (where R' 

is the c — B relative coordinate), and represents a 
scattering state for the x + T system. By making use 
of the Gell-Mann-Goldberger two-potential formula, the 
transition amplitude 0 can be rewritten as: 

Tpost = (XcrVir Ii4:c14'j- ^), (5) 

where xir is distorted wave generated by the poten¬ 
tial Uc- The above matrix element is dominated by small 
X — c separations, where jg j^g 

Although very appealing from the formal point of view, 
expression 0 is hard to implement in practice. The main 
reason is that this expression involves a six-dimensional 
integral, in which both the initial and final wavefunctions 
are highly oscillatory. Furthermore, post form represen¬ 
tations offer poor convergence since both the scattering 
waves for a; -I- c and the potential 14c are expressed in 
the same coordinate and consequently there maybe no 
natural cutoff for the integral 0 Hi- 
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Figure 2: : (a) Breakup couplings for a two-body projectile 
(c+x) impinging on a target T and (b) corresponding transfer 
to the continuum couplings. 

In order to make the calculation more feasible, Shyam 
and collaborators (see, for instance, [^) have 

developed an approach based in the amplitude 0, in 
which the exact wavefunction is replaced by its elastic 
component, i.e., 

vI/(+)«x^+HRAT)0o(r.c), (6) 

where (paif^xc) is the projectile ground state and is 
a distorted wave generated with a optical potential, typi¬ 
cally adjusted to reproduce the elastic scattering data. 
Note that this approximation neglects breakup events 
that proceed via projectile excitation. Note also that, 
even after the replacement of the exact wavefunction by 
the elastic component, Eq. 0 still involves a six dimen¬ 
sional integral. A significant simplification of the prob¬ 
lem can be achieved by using the local momentum ap¬ 
proximation [^, which leads to a factorization of 
the amplitude in a product of two terms, each involving 
a three-dimensional integral. 

The difficulties outlined above can be partially avoided 
using the prior representation of the transition ampli¬ 
tude, 

Tprior = ('kj ''iVxT + Uc — UAxlfpOXo), (7) 

where the final state is the exact three-body wavefunction 
with incoming boundary conditions. This wavefunction 
is typically expanded in terms of the x + T continuum 
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states, formally similar to Eq. but now in the Jacobi 
set (2) of Fig.n 

N2 

= ( 8 ) 

0=0 

Again, the functions (fi/s represent the set of bin wave- 
functions umiiii, constructed by superposition of 
pure scattering waves. Note that Eq. m goes beyond 
the DWBA, because couplings between final states are 
explicitly considered in the wavefunction 4'} . 

Therefore, in the standard CDCC method the three- 
body continuum is described in terms of the projectile 
two-body {x — c) states, while in the amplitude 0 this 
continuum is expanded using the fragment-target states 
x — T. While the CDCC method treats the breakup pro¬ 
cess as inelastic excitations to the projectile continuum, 
expressions m and 0 emphasize a rather different pic¬ 
ture, in which three-body breakup is formally treated as 
transfer of one of the fragments {x in our case) to the 
target continuum. This is schematically illustrated in 
Fig.EIb). At this stage, it is worth to stress that, in the 
way here presented, the CDCC and the the TR* meth¬ 
ods are solutions of the same three three-body model 
Hamiltonian, given by Eq. 0 and, therefore, three-body 
observables obtained with these two approaches should 
be the same. 

However, in practice there are several factors that may 
destroy this equivalence. First, due to computational lim¬ 
itations, one can not include an arbitrarily large number 
of continuum states. Secondly, there are ambiguities as¬ 
sociated with the choice of the interactions involved in 
both schemes. In the BU approach, one usually has two 
complex potentials, namely Ux and 14, and a real in¬ 
teraction, Vxc- By contrast, in the amplitude 0 the 
wavefunction 4'), ^ is typically obtained with the complex 
potentials Uc and Uxc- The choice of the potential I4 t 
deserves special care. For inclusive processes, in which 
the fragment x is allowed to interact in any possible way 
with the target, VxT would be a complicated many body 
operator, which can induce excitations in both x and T. 
However, for a comparison with CDCC, in which only 
the elastic breakup component is calculated, this opera¬ 
tor is better represented by an effective complex optical 
potential miHi and hence, according to our previous 
notation, f4T = Ux- This is actually the choice made in 
current semiclassical applications of the transfer to the 
continuum method 0. 

Another ambiguity is related to the interaction that 
should be used for Uat in the amplitude 0 . Note that, 
if the exact expression is used for the three-body wave- 
function ^ the matrix element is independent of the 
choice of the potential Uat- This result does not hold 
when ^ is replaced by an approximated wavefunction. 
Following the standard DWBA choice, one could use the 
optical potential that reproduces the elastic scattering. 
Another possible choice, is the so called cluster-folding 


potential, given by the sum of the fragments-target in¬ 
teractions folded with the ground state of the projectile: 
{(j^olUx + Uc\(j)o)- In our calculations we have explored 
both choices. 

The main purpose of this work is to test to what ex¬ 
tent the equivalence between the BU and (prior form) 
TR* is satisfied, at least in an approximate way, in ac¬ 
tual calculations. To this end, we have performed nu¬ 
merical calculations for two different systems using both 
approaches, and compared several reaction observables. 
At high scattering energies, around 100 MeV per nucleon 
and above, the TR* method, as presented here, becomes 
numerically very demanding, and the problem is better 
solved by using further approximations, such as the use 
of classical trajectories. Since it is our purpose to com¬ 
pare the full quantum mechanical CDCC and TR* ex¬ 
pressions, we confine ourselves to reactions at low and 
medium energies. 


III. CALCULATIONS 

A. p-|-^^Be case 

We first consider the reaction of 38.5 MeV per nucleon 
^^Be breaking up on protons. The elastic and trans¬ 
fer channels were measured in GANIL i HI but no 
breakup data was recorded. According to the discus¬ 
sion in the previous section, the ^^Be breakup reaction 
can be thought of as the direct breakup (BU) ^^Be-bp ^ 
(^°Be-|-n) +p or transfer of the neutron to the continuum 
of the deuteron (TR*) ^^Be-bp ^^°Be-b(n + p). 

The n — p interaction was taken from [l^, whereas 
the nuclear interaction for p—^°Be was extracted from a 
fit to the elastic data j^. The Coulomb potential for 
p—^°Be was also included so Coulomb breakup is also in¬ 
cluded in our calculations, although it was shown to be 
very small. The binding potential and the potential gen¬ 
erating the continuum waves for n—^°Be was the same 
as in Q, but without the spin-orbit term. These poten¬ 
tials are listed in Tabled The BU calculations required 
partial waves up to Imax = 4 and energies up to £max=30 
MeV for the relative motion of the n—^°Be system. The 
bin wavefunctions for the CDCC couplings were calcu¬ 
lated up to Rbiri=60 fm. An Umax = 25 was necessary 
for the ^^Be—p distorted waves. As to the TR* calcula¬ 
tion, the same parameters were sufficient for convergence 
but they are computationally more lengthy ji^. All the 
TR* calculations here presented use as incoming optical 
potential the folding of the p — n and p-^°Be interactions 
with the ground state wavefunction of the ^^Be nucleus. 
We also did calculations using a Woods-Saxon shape with 
the same parameters as for the p-^°Be potential. Results 
obtained with this potential are very similar to those of 
the cluster-folding and, hence, will not be shown in the 
graphs. Both the BU and TR* calculations were per¬ 
formed with the computer code FRESCO 

In Fig. 0 we present the differential breakup cross sec- 
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Table I: Optical model parameters used in this work. Except 
for the p + n case, all potentials are parameterized using the 
usual Woods-Saxon form, with a real volume part and vol¬ 
ume (Wv) or surface (Wd) imaginary part. Reduced radii are 
related to physical radii by i? = . 


System 

Vo 

(MeV) 

ro 

(fm) 

0-0 

(fm) 

1F„ 

(MeV) 

Wd 

(MeV) 

Ti 

(fm) 

0,1 

(fm) 

p+^°Be 

51.2 

1.114 

0.57 

19.5 

0 

1.114 

0.50 

p-f n“ 

72.15 

1.484 

- 

- 

- 

- 

- 

®B+®®Ni 

130 

1.050 

0.65 

92 

0 

1.123 

0.997 

■^Be+^^Ni 

100 

1.050 

0.65 

30.6 

0 

1.123 

0.80 

p+®«Ni 

54.512 

1.17 

0.75 

0 

11.836 

1.260 

0.58 

p+'^Be ’’ 

44.675 

1.25 

0.52 

- 

- 

- 

- 


“Gaussian geometry: V(r) = Vb exp[(r/ro)^]. 

^In the TR* case, this potential includes also an imaginary part 
with the same geometry as the real part. 



Figure 3: (Color online) Breakup energy distribution for ^^Be 
on protons at 38.5 MeV/u: comparison of transfer to the con¬ 
tinuum (bottom figure) with the direct breakup approach (up¬ 
per figure). 


tion calculated within the BU and TR* methods, as a 
function of the excitation energy of the ^°Be-n and p- 
n systems, respectively. It can be seen that, in the BU 
case, most of the strength is below e^iC^^Be) « 5 MeV 
whereas in the TR* calculation the strength is largely 
concentrated around Sxid) « 20 MeV. The total inte- 



Figure 4: Total angular momentum distribution for the 
breakup of ^^Be on protons at 38.5 MeV/u: comparison of 
transfer to the continuum with the direct breakup approach. 


grated cross section for the two processes are abu = 125 
mb and atr = 140 mb. 

In order to establish a meaningful comparison between 
the two approaches, one has to compare the same quanti¬ 
ties. For this purpose, in Fig. 0 we compare the contribu¬ 
tion of each total angular momentum J (resulting from 
the vector coupling of projectile, target and their relative 
motion angular momentum) to the breakup cross section. 
It can be seen that both distributions are similar for small 
values of J. Also, we find that for the two cases the distri¬ 
bution peaks around J = 5 which means that most of the 
breakup cross section occurs at distances 6 4 fm. The 

similitude between both distributions supports the idea 
that breakup, calculated as excitation of the projectile 
to its continuum spectrum, or by transfer to the contin¬ 
uum states of the target, do describe the same physical 
process. However, for J > 5 the TR* clearly exceeds the 
BU cross section which, as we will show below, results on 
different predictions for measurable physical observables. 

In actual breakup experiments, the data commonly 
recorded are the angular and/or energy distributions of 
the emerging fragments. Therefore, it is instructive to 
compare the predictions of both approaches for these ob¬ 
servables. In Fig. 13 we represent the calculated breakup 
cross section distribution of the outgoing ^°Be fragments 
as a function of its kinetic energy, measured in the overall 
c.m. of the three-body system. In both methods, these 
distributions are obtained by integration of a triple dif¬ 
ferential cross section with respect to the angular vari¬ 
ables. In the case of the TR* approach, this procedure is 
straightforward, since expression 0 is referred already to 
the scattering angle and energy of the ^°Be fragments. In 
the case of the BU approach, the differential cross section 
is naturally expressed in terms of the scattering angle of 
the composite x + c and the relative energy between these 
two fragments. In order to obtain the differential cross 
section with respect to any of the fragments, one has to 
apply t o ap propriate kinematic transformation, as done 
in Ref. )28l |. 

These energy distributions show only a qualitative 





























6 



Figure 5: Energy distribution in the overall c.m. for the 
coming from the breakup of on protons at 38.5 
MeV/u: comparison of transfer to the continuum with the 
direct breakup approach. 


agreement between the two methods. In both cases, the 
energy of the ^°Be fragments goes from zero to about 6 
MeV, with a maximum around 3 MeV. However, although 
the same energy region of space is being included, the 
two models do not produce identical shapes. In Fig. [3 we 
also show the convergence rates for both TR* and BU. 
The labels s,p,d, etc refer to the relative partial waves 
included in the corresponding calculation. For instance, 
the solid line in the BU calculation includes all n—^°Be 
partial waves up to ^ = 4 whereas the dashed line includes 
only partial waves up to ^ = 2. 

Disagreements become more severe for the angular dis¬ 
tributions. These can be seen in Fig. |H1 Note that the 
TR* calculation (dot-dashed curve in Fig. EJ exhibits a 
pronounced decrease of the cross section as a function of 
angle. In addition, its forward angle cross section is an 
order of magnitude larger than the BU calculation (solid 
line) and an order of magnitude lower for backward an¬ 
gles. Part of the reason for the disagreement can be un¬ 
derstood excluding the d-wave resonance in the n—^°Be 
system (dashed curve). This wave has a very strong con¬ 
tribution for backward angles and a d-wave resonance in 
^^Be will be very hard to model in terms of the deuteron 
continuum. However, the discrepancy remains at forward 
angles: including only the s-wave of the deuteron in the 
BU calculation, the resulting cross section is an order of 
magnitude smaller than the TR* cross section. Detailed 
data for this reaction would be very useful. 

B. case 

A breakup reaction for which more detailed data exist 
is that for ®B^’^Be-|-p on ^®Ni at 25.6 MeV. Calculations 
using the standard CDCC ®B-|-®®Ni ^ (^Be-|-p)-l-^®Ni 
ha ve p rovided very good agreement with experiment 
0, [ 2 ^. Again, one can think of the alternative path 
to breakup, as transfer to the continuum of the ^®Cu nu¬ 
cleus (TR*) ®B-|-®®Ni ^^Be-|-(p-|-®®Ni). All interactions 



Figure 6: Angular distribution in the c.m. for the ^°Be com¬ 
ing from the breakup of ^’^Be on protons at 38.5 MeV/u: com¬ 
parison of transfer to the continuum with the direct breakup 
approach. 


for both BU and TR* are the same as those in [^, al¬ 
though for the p—Ni only the real part was included in 
the TR* calculation. 

The BU calculations required partial waves up to 
^max = 4 and energies up to Smax = 8 MeV for the rela¬ 
tive motion of the p-^Be system. The bin wavefunctions 
for the CDCC couplings were calculated up to Rbin = 60 
fm and the coupled channel equations were solved with 
Rmax = 500 fm. An Lmax = 1000 was necessary for 
the ®B—®®Ni distorted waves ^3|. Note that in this 
case the BU calculation is a good test reference, as it 
agrees ver y w ell with both energy and angular distribu¬ 
tion data 1^, at least within the kinematic conditions 
of the referred experiment. As to the TR* calculation 
we used partial waves up to lmax = 17 and energies 
up to Emax = 10 MeV for the relative motion p-|-®®Ni. 
To reduce the computational requirements, for Z/ > 6, 
continuum-continuum couplings were included only be¬ 
tween bins with the same If. The bin wavefunctions were 
calculated up to Rbin = 120 fm, and an Lmax = 120 
was necessary for the ®B—^®Ni distorted waves. How¬ 
ever, these results are not yet converged. The large re¬ 
quired widths for the non-local transfer couplings make 
the calculations extremely heavy. 

We next compare the same quantities as in the p+^^Be 
case. Unlike the previous test example, where Vxc = Vpn, 
here the I4c interaction (p+^Be), as extracted from the 
elastic data, is expected to contain an imaginary part. 
In our TR* calculations, we probe several possibilities 
for the imaginary part, keeping the same geometry as 
the real part and using different choices for the depth. 
For the incoming channel optical potential we used two 
different potentials. The first one, denoted OMl, is the 
sum of the p+^Be and ^Be+^®Ni interactions folded with 
the bound state wavefunction of the ®B nucleus. The sec¬ 
ond one, consisted on a parametrization with two Woods- 
Saxon terms, real and imaginary, with parameters ob¬ 
tained by fitting the elastic angular distribution, as pre- 
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Figure 7: Total angular momentum distribution for the 
breakup of ®Be on ®®Ni at 25.6 MeV. The solid line is the 
BU calculation, whereas the remaining lines correspond to 
the TR* calculations with two different sets of partial waves 
lf,as indicated by the labels (see text for details). 



Figure 8: (Color online) Angular distribution in the labora¬ 
tory frame for the ^Be fragments coming from the breakup of 
®B on ®®Ni at 25.6 MeV within the direct breakup and trans¬ 
fer to the continuum approaches. Experimental data are from 


dieted by the CDCC calculation. In this case we used, as 
starting parameters for the fitting routine, those for the 
^Be+®®Ni interaction. The parameters for this potential 
resulting from the fit, denoted OM2, are listed in Table 

III 

In Fig. [3 we show the total breakup calculated within 
the BU and TR* schemes, for each value of the total an¬ 
gular momentum, J. In the TR* case, two different cal¬ 
culations are presented. In both cases, the potential OM2 
was used for the elastic channel, and the value Wd=3 was 
used for the imaginary depth of the p+^Be interaction. 
The thin solid line in this figure represents the TR* cal¬ 
culation performed in the subspace If = 0 — 6. This 
calculation exhibits clear differences from the BU (thick 
solid line): the lower values of J are clearly overesti¬ 
mated, whereas for the large values of J, the distribution 
falls too fast as compared to the BU. The second TR* 
calculation here presented (dotted-dashed line) uses the 
same parameters as before, but includes also the partial 
waves Z/ = 7 — 17 for the proton-®®Ni continuum. This 
calculation improves the agreement for large J. However, 
it also adds an extra contribution on the lower values of 
J which appears to deteriorate the agreement with the 
BU results. 

As noted in the previous section, it has been argued 
that in the calculation of elastic breakup, one should use 
an absorptive potential for the I4 t operator in Eq. 0 . 
This might explain part of the discrepancies found here 
between both approaches. Unfortunately, the present 
version of the code FRESCO does not allow this interac¬ 
tion to be complex. Although at present we cannot use 
an imaginary term in the 14t operator, we are consider¬ 
ing modifications of the code to enable this feature. R 
is fortunate that, in the p+^^Be reaction, Vxt — Vpn 
which, within the energy range of our analysis, is well 
represented by a real potential. This might explain the 
better agreement obtained for the absolute value of the 
breakup cross section, as well a for the J-distribution, as 


compared to the ®B+^®Ni case. 

Although the representation of the Vxt operator by a 
real operator has undesirable consequences for the pur¬ 
pose of the present work, one may speculate about the 
physical meaning behind this choice. As noted above, 
choosing I4 t as the phenomenological optical potential 
in the CDCC approach implies that the model space con¬ 
siders only breakup where the target is left in its ground 
state. Conversely, one may argue that, by choosing this 
potential as real, as we do here in the TR* case, there is 
the possibility of including other inelastic breakup events, 
and even transfer to bound states of the target. In other 
words, the TR* method with a real Vxt may include con¬ 
tributions from the so called stripping breakup. These 
aspects will be explored in future work. 

For the ^Be angular distribution, we first study the 
convergence of the calculations with respect to the size 
of the model space. In the case of the BU this has been 
discussed in detail in Ref. [i^, and so we will just quote 
the results from this reference. In this section, we concen¬ 
trate only on the convergence of this observable within 
the TR* scheme. For definiteness, these calculations were 
performed with the choice Wd = 3 MeV for the imagi¬ 
nary part of the p+^Be potential. The dependence on 
this potential will be analyzed below. The convergence 
of the calculation with respect to the number of partial 
waves for the relative motion of the x — T pair is depicted 
in Fig.lSI For comparison, the BU calculation (thick solid 
line) and the experimental data points have been also 
included. The thin solid line is the calculation where 
the partial waves If = 0 — 6 are included and coupled 
among them to all orders. The thick dotted-dashed line 
is the sum of this calculation and the separated differen¬ 
tial cross sections for the partial waves If = 7 — 17. It 
becomes apparent that the contribution of these partial 
waves is very important to describe the strong Coulomb 
peak at small scattering angles. Both calculations are in 
good agreement with the data. Given the large error bars 
and restricted angular range of the data it is not possi- 
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Figure 9: Angular distribution in the c.m. for the ^Be 
fragment coming from the breakup of ®B on ®®Ni at 25.6 
MeV: comparison of transfer to the continuum with the direct 
breakup approach for two different choices of the incoming op¬ 
tical potential for the TR* calculation (see text for details). 


ble to make strong conclusions on which method is more 
suitable in this particular situation. Roughly speaking, 
it seems that the TR* describes better the larger angles, 
while the BU is more suitable to describe the smaller an¬ 
gles. In the remaining discussion, all our comparisons 
with the BU calculations will be performed with the full 
set of partial waves {If = 0 — 17). 

Next, we study the dependence of the ^Be angular dis¬ 
tribution on the choice of the incoming channel optical 
potential. This is illustrated in Fig. El In this figure, the 
dashed and dotted-dashed lines correspond, respectively, 
to the TR* calculation with the cluster-folded potential 
(OMl) and the phenomenological optical potential ob¬ 
tained from a fit of the CDCC elastic angular distribu¬ 
tion (OM2). Again, we have fixed the imaginary depth 
of the p+^Be interaction to Wd = 3 MeV. For compar¬ 
ison purposes, the BU calculation (thick solid line) has 
been also included. One sees that the choice of the elas¬ 
tic channel optical potential has indeed an effect on the 
predicted ^Be cross sections. However, given the uncer¬ 
tainties of these calculations, differences do not seem very 
dramatic, and in both cases a fairly good agreement is ob¬ 
tained with the BU calculation irrespective of the choice 
of this potential. 

In Fig. El we compare the standard CDCC calculation 
with TR* calculations performed with the optical poten¬ 
tial OM2 for the incoming channel, and different choices 
of the imaginary depth of the p+^Be potential. Thick 
lines correspond to the full calculation {Imax = 17). At 
backward angles all TR* calculations look very similar, 
indicating a fast convergence with respect to the number 
of partial waves and a weak dependence on the choice of 
the p+^Be potential at these angles. The effect of the 
imaginary part seems to be crucial at intermediate an¬ 
gles, where one observes a progressive suppression of the 
cross section with increasing absorption. The TR* with 
real p+^Be interaction clearly overestimates the BU re¬ 
sult at intermediate angles. Interestingly, the forward 
angular region is only weakly affected by this absorptive 



Figure 10: (Color online) Angular distribution for the ^Be 
fragments coming from the breakup of ®B on ®®Ni at 25.6 
MeV: comparison of transfer to the continuum with the direct 
breakup approach for different values of the imaginary depth 
of the p+^Be potential. 


term. As we verified in our calculations, this is a conse¬ 
quence of the fact that this potential has little effect on 
the higher partial waves. The best agreement with the 
BU calculation is obtained when the imaginary part of 
the p+^Be potential for TR* is Wd = 6 MeV. 

The reasonably good agreement between BU calcula¬ 
tion and the TR* calculations, performed in the aug¬ 
mented model space {Imax = 17) and with a complex 
proton+'^Be interaction, leads us again to the conclusion 
that projectile breakup and transfer to the target con¬ 
tinuum populate, to a large extent, the same three-body 
continuum. We interpret the discrepancy at small scat¬ 
tering angles as lack of convergence of our TR* calcula¬ 
tions, and the ambiguities associated to the potentials. 

In Fig. El we plot the energy distribution for the de¬ 
tected ^Be fragments after breakup, for the BU (solid 
line) and the TR* approach with Imax = 17, and Wd = 3 



Figure 11: Energy distribution in the c.m. for the ^Be coming 
from the breakup of ®B on ®®Ni at 25.6 MeV: comparison of 
transfer to the continuum with the direct breakup approach. 
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Figure 12: Breakup distribution as a function of the relative 
motion ®®Ni-p for the transfer to the continuum of the proton 
from ®B to ®®Ni at 25.6 MeV. Results up to // =6 are fully 
coupled, but for 1/ >6 they are calculated separately (see 
text for details). 

MeV for p+'^Be (dashed line). It can be seen that both 
methods give similar distributions. In particular, it is 
noticeable that both methods predict a maximum of the 
energy distribution at about the same ^Be energy. The 
TR* calculation gives however a larger breakup cross sec¬ 
tion. 

To have further insight into the convergence of the TR* 
calculation with respect to the size of the model space 
we have plotted in Fig. El the distribution of the TR* 
cross section as a function oi If, i.e., the final angular 
momentum between the proton and the target. On one 
hand, it is clear that TR* requires far more partial waves 
that the BU calculation. On the other hand, the small 
contribution for > 15, does provide some confidence in 
the results presented here. The breakup of ®B on ®®Ni at 
25.8 MeV is a good example where the BU configuration 
seems to work better than the TR* configuration. 

IV. DISCUSSION 

The qualitative agreement between the calculations 
performed in the BU and TC representations clearly indi¬ 
cates that both basis describe to a large extent the same 
three-body continuum. However, the analysis of the pre¬ 
ceding section also shows that, in order to achieve con¬ 
vergence of the observables, the number of basis states 
required in both representations can be very different. 
From the practical point of view, it will be desirable in 
general to choose the representation that requires less 
number of states. 

For example, our analysis of the ®B breakup reaction 
clearly supports the choice of the Jacobi coordinates (1) 
in FigO] On the other side, our previous study on the 
reaction ®Li+^°®Pb was better performed using the 
coordinate set (2). It is clear that, in general, the most 
suitable choice will depend on the specific reaction. In¬ 
spired by the work of Merkuriev on three-body bound 
states, we have searched for a criterion that can select be¬ 


tween the two representations. Unfortunately the asymp¬ 
totic behaviour of the three body continuum is very dif¬ 
ferent from the exponential decay of bound states, and 
the final behaviour of these expansions is not as trans¬ 
parent. Therefore, we will present only qualitative argu¬ 
ments to evaluate the relative importance of the different 
configurations. 

For practical purposes one always opts for the calcula¬ 
tion that requires the minimum number of partial waves 
in the x—c or x—T systems, for BU and TR* respectively. 
For the ^^Be the difference in Imax for BU and TR* was 
not noticeable. For the ®B example, Imax = 4 was suf¬ 
ficient for BU whereas Imax = 17 was still not enough 
for the TR*. In addition to the angular momentum, it 
seems clear that if a representation such as Eq. 0 is 
valid, then the average energy {ei) associated with the 
relative coordinate ri should be much less than the to¬ 
tal energy in the centre of mass frame Ecm. Equally, 
if Eq. (|Hl is to be used, the average energy (£ 2 ) associ¬ 
ated with coordinate r 2 should be small compared to the 
total energy of the exit partition in the centre of mass 
frame Ecm- We have computed the average energy be¬ 
tween the fragments in the continuum, weighting it by 
the cross section. For the ®B example above, we obtain 
(ei) = 1.85 MeV with E^cm = 22.7 MeV for the BU cal¬ 
culation and (£ 2 ) = 7.84 MeV with E^m = 26.0 MeV for 
the TR* calculation. In the first case (ei)/Ecm = 0.08 
whereas in the latter ( 62 )/Ecm = 0.30. As to the ^^Be 
example, the difference is also pronounced: for BU we 
obtain (ei)/ e/I/ = 0.15 and {e 2 )/Eim = 0.52, imply¬ 
ing again that the transfer to the continuum approach is 
not the best. In the breakup of ®Li [^, the TR* ap¬ 
proach was used successfully. We compute the average 
energy between the fragments in the outgoing channel 
and obtain ( 89 )/E/m = 0.05, validating the previous TR* 
calculations [^. This same criterion shows a red card 
to the preliminary calculations on ®He [^, as in that 
case we have /e’^/E/m = 0.43. In conclusion, the con¬ 
dition <^C 1 should be satisfied whenever only 

one i Jacobi set is taken into account in the reaction 
formalism. This is for instance the case in Coulomb dis¬ 
sociation, where the long-range and smooth behaviour of 
the Coulomb potential makes that the reaction mecha¬ 
nism populates mainly low-energy states of the projec¬ 
tile. As a matter of fact, is was shown in 0 that the 
peak in the breakup angular distribution at 0c.m. ~ 15° 
is mainly due to Coulomb excitation. In our calculations, 
this peak is well reproduced by the CDCC calculation, 
while it requires many partial waves of the proton-target 
system in the TR* method. Furthermore, at high en¬ 
ergies, the Coulomb dissociation cross section becomes 
approximately proportional to the B(E1) strength of the 
projectile. Hence, the BU approach provides in this case 
a more transparent and useful picture of the reaction pro¬ 
cess. On the other side, in situations where the removed 
particle has a high probability to be left with a small 
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relative energy with respect to the target, the TR* may 
result more convenient. Nevertheless, this does not seem 
to be the case of the reactions studied in this work. 

Even in calculations where the TR* converges quickly 
with the number of partial waves, we found it to be 
computationally very demanding, as it involves non-local 
couplings. In practice, these calculations could be signif¬ 
icantly speeded up, using different techniques which are 
now of common use: local momentum and adiabatic ap¬ 
proximations, etc. It is however beyond the scope of this 
work to explore how to make the TR* numerically more 
feasible. 


V. CONCLUSIONS 

Given the importance of the reaction model in the un¬ 
derstanding of the fundamental nuclear structure on the 
drip lines, we compare two alternative schemes to calcu¬ 
late breakup observables for the reaction A{= c + v)+T, 
within the same three-body Hamiltonian. Each one of 
these methods uses a description of the three-body con¬ 
tinuum in terms of one of the possible sets of Jacobi co¬ 
ordinates. In the CDCC approach, the three-body con¬ 
tinuum is described in terms of the c — v states. On the 
contrary, the transfer to the continuum (TR*) approach 
expands the continuum in terms oi v — T states. Since 
both sets of states form a complete basis, reaction ob¬ 
servables could be in principle calculated using either of 
these two basis. We show that in both cases predictions 
by these two schemes are in semiquantitative agreement. 
This result clearly shows that, provided that enough basis 


states are included, both representations describe essen¬ 
tially the same three-body continuum. In the 
reaction, both calculations are consistent with existing 
experimental data. From our analysis, it is clear that 
the truncated model space is not always identical. In 
particular, in this case we found that the TR* approach 
requires a significantly larger number of partial waves for 
the proton-cluster relative motion, thus making the cal¬ 
culation numerically more demanding. We also find that 
part of the disagreement between the two methods is due 
to the ambiguities associated with the choice of the effec¬ 
tive interactions involved in both methods. In addition, 
we have proposed a simple criterion based on the aver¬ 
age relative excitation energy to select between the two 
approaches. 

More detailed studies on the absorption part of the op¬ 
tical potential and applications to other reactions will be 
presented elsewhere. Ultimately, we would like to com¬ 
pare these results with exact Faddeev calculations. Work 
along these line is being initiated. 
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